Precision Measurements of Collective Oscillations in the BEC-BCS Crossover 
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We report on precision measurements of the frequency of the radial compression mode in a strongly 
interacting, optically trapped Fermi gas of 6 Li atoms. Our results allow for a test of theoretical 
predictions for the equation of state in the BEC-BCS crossover. We confirm recent quantum Monte- 
Carlo results and rule out simple mean-field BCS theory. Our results show the long-sought beyond- 
mean- field effects in the strongly interacting BEC regime. 
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Ultracold, strongly interacting Fermi gases [l|, 0, 0, 0, 
0, 0, 0, 0, U 03 111 H, El have attracted consider- 
able attention over the past few years, serving as unique 
model systems to create, control, and investigate novel 
states of quantum matter. Experimentally, the availabil- 
ity of such systems has opened up exciting possibilities 
to study many-body quantum phenomena like molecular 
Bose-Einstein condensation (BEC) [1 and the crossover 
from BEC to a Bardeen-Cooper-Schrieffer (BCS) type 
superfluid 0, i, 0, 0, 0, 0, M, [H [U, El • These experi- 
ments may also lead to a better understanding of strongly 
interacting quantum systems in different areas of physics, 
ranging from high-T c superconductors to neutron stars 
and the quark-gluon plasma. 

A degenerate two-component Fermi gas undergoes the 
BEC-BCS crossover [lj] when the s-wave scattering 
length a is varied from positive to negative values across 
a scattering resonance. In the crossover region, where a 
is comparable with or larger than the interparticle spac- 
ing, the equation of state is governed by many-body ef- 
fects. Understanding the equation of state is a funda- 
mentally important challenge and constitutes a difficult 
task for many-body quantum theories, even in the zero- 
temperature limit. Mean-field BCS theory [l4| provides 
a reasonable interpolation between the well-understood 
limits. More sophisticated crossover approaches [l5| yield 
quantitatively different results in certain regimes, none of 
them however providing a complete description of the 
problem. The most advanced theoretical results were 
obtained by numerical calculations based on a quantum 
Monte-Carlo (QMC) approach [ll]. 

On the BEC side of the crossover, there is an in- 
teresting competition in the equation of state between 
the strong interactions in a Bose gas and the onset of 
fermionic behavior. For a strongly interacting Bose gas, 
one can expect quantum depletion to increase the average 
energy per particle. To lowest order, this beyond-mean- 
field effect leads a correction to the equation of state pre- 
dicted by Lee, Huang, and Yang (LHY) almost 50 years 
ago Beyond mean-field effects are expected to re- 

duce the compressibility of a strongly interacting Bose 
gas as compared to the weakly interacting case. How- 



ever, when approaching the resonance, fermionic behav- 
ior emerges and the system loses its purely bosonic char- 
acter, which increases the compressibility of the strongly 
interacting gas. Mean-field BCS theory does not con- 
tain beyond-mean field effects and the LHY correction is 
absent there. However, the QMC results predict beyond- 
mean-field effects to be visible on the BEC-side of the 
crossover [l6| . 

In this Letter, we report on precision measurements 
of the radial compression mode in an optically trapped, 
strongly interacting Fermi gas of 6 Li atoms. The mode 
serves as a sensitive probe for the compressibility and 
thus the equation of state of a superfluid gas in the BEC- 
BCS crossover [TH, [II . We reach a precision level that 
allows us to distinguish between the predictions result- 
ing from mean field BCS theory and QMC calculations. 
Previous experiments on collective modes, performed at 
Duke University [j| [ll| and at Innsbruck University d , 
showed frequency changes in the BEC-BCS crossover in 
both the slow axial mode and the fast radial compression 
mode of a cigar-shaped sample. The accuracy, however, 
was insufficient for a conclusive test of the different many- 
body theories in the strongly interacting regime. 

We prepare a strongly interacting, degenerate gas of 
6 Li atoms in the lowest two internal states as described 
in our previous publications [3, 0, [l(| • The broad Fesh- 
bach resonance centered at a magnetic field of B = 834 G 
facilitates precise tuning of the scattering length a [11 . 
Forced evaporative cooling is performed in a 1030-nm 
near-infrared laser beam focussed to a waist of 54/im 
at 764 G. This results in a deeply degenerate cloud of 
iV = 2.0(5) x 10 5 atoms. By adiabatically increasing 
the trap laser power after cooling, the sample is recom- 
pressed to achieve nearly harmonic confinement. In the 
axial direction the gas is magnetically confined in the cur- 
vature of the field used for Feshbach tuning with an axial 
trap frequency of w z /2ir = 22.4 Hz at 834 G. The experi- 
ments reported here are performed at two different final 
values of the laser power of the recompressed trap. At 
135 mW (540mW), the trap is 1.8^K (7.3/xK) deep and 
the radial trap frequency is uo r « 27r x 290 Hz (590 Hz). 
The Fermi energy of a non-interacting cloud is calcu- 
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FIG. 1: Radial sloshing observed at a trap power of 540 mW 
and B = 735 G (1/kpa = 1.55). The two-dimensional center- 
of-mass motion is represented in a coordinate system (x',y') 
rotated by 45° with respect to the principal axes of the trap. 
The beat signal between the two sloshing eigenmodes demon- 
strates the ellipticity of the trap with the two eigenfrequencies 
w x /27T = 570 Hz and uj y /2n = 608 Hz (ellipticity e = 0.066). 



lated to E F = h 2 k 2 F /2m = h(3w*w z Ny/ 3 = k B x 500 nK 
(800 nK); here m is the mass of an atom and k B is Boltz- 
mann's constant. 

Since our first measurements on collective excitation 
modes @, we have upgraded our apparatus with a two- 
dimensional acousto-optical deflection system for the 
trapping beam and a new imaging system along the trap- 
ping beam axis. These two improvements provide us with 
full access to manipulate and observe the radial motion. 

The trap beam profile is somewhat elliptic because of 
imperfections and aberrations in the optical set up. To 
simultaneously excite the two eigenmodes of the radial 
sloshing motion, we initially displace the trapped sample 
into a direction between the horizontal and vertical prin- 
cipal axes of the radial potential. After a variable hold 
time, during which the cloud oscillates freely we turn 
off the optical trap. After a time of flight of typically 
4 ms we take an absorption image of the released cloud. 
The center-of-mass position of the cloud then reflects its 
momentum at the instant of release. The experimental 
results in Fig. [1] demonstrate the sloshing with a beat 
between the two eigenmodes. A careful analysis of such 
data [2l| allows us to determine the eigenfrequencies u> x 
(horizontal sloshing) and w y (vertical sloshing) to within 
a relative uncertainty of typically 2 x 10~ 3 . We finally 
derive the mean sloshing frequency u>j_ — ^/bj x uj y and the 
ellipticity parameter e = {oJ y — lu x )/lo±. 

To excite the radial compression oscillation we reduce 
the trap light power for a short time interval of ~100/xs, 
inducing an oscillation with a relative amplitude of typi- 
cally 10%. After a variable hold time the cloud is released 
from the trap. From fits of two-dimensional Thomas- 
Fermi profiles to images of the expanding cloud taken 
4ms after release, we determine the mean cloud radius. 
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FIG. 2: Radial compression oscillation observed for the same 
conditions as the sloshing mode data in Fig. [1] The radial 
width is determined by averaging the horizontal and vertical 
Thomas-Fermi radii after expansion. Here we obtain u> c /2tt = 
1185 Hz. 

A typical set of measurements is shown in Fig. [2j A fit 
of a damped harmonic oscillation to such data yields the 
frequency lo c and damping rate 7 of the radial compres- 
sion mode. 

Our experiments are performed close to the limit of an 
elongated harmonic trap potential with cylindrical sym- 
metry. This elementary case is of great general relevance 
for many quantum gas experiments in optical and mag- 
netic traps (see, e.g., [22| ), and collective excitations are 
conveniently normalized to the trap frequency u> r (TR[l9j. 
The compression mode frequency can then be written as 
w c = f c ijJ r , where f c is a dimensionless function of the 
interaction parameter 1/fcpa and is related to an effec- 
tive polytropic index T [H, [l9| of the equation of state 

by ^ c 2 = 2(r + iK 2 . 

In order to compare our experimental results with the- 
ory we consider the quantity f c , i.e. the normalized com- 
pression mode frequency of the ideal, cylindrically sym- 
metric, elongated trap. We find, that for our experimen- 
tal conditions, f c is approximated by the ratio u c /luj_ of 
the measured compression mode (lu c ) and mean slosh- 
ing mode frequencies to better than one percent. 
On the desired accuracy level of 10~ 3 , however, two 
small effects have to be taken into account: the resid- 
ual trap ellipticity and the anharmonicity of the radial 
potential in combination with the spatial extension of 
the trapped sample. We thus introduce two small cor- 
responding corrections, expressing f c in the form f c — 
(1 — ne 2 + ba) uj c /lu±. 

For the ellipticity correction ne 2 , a straightforward so- 
lution of the hydrodynamic eigenfrequency equation [2lT | 
yields k = (2 + T)/4T, where T can be approximated by 
r = (u) c /l>±) 2 /2 — 1. For the anharmonicity correction, 
the parameter a — ^muj 2 _r 2 nls /Uo relates the potential 
energy associated with the root-mean-square radius r rms 
[23l ] of the trapped cloud to the trap depth. The coeffi- 
cient b results from the differential anharmonicity shifts 
in the compression and sloshing modes and can be calcu- 
lated according to [HI, 0, [H] . We obtain 6 = 0.167 
and 0.280 in the limits of BEC and unitarity, respectively. 

Our measurements on the sloshing and compression 
modes are summarized in Table U including the two 
small corrections. For the data in the strongly interacting 



3 



TABLE I: Experimental data on radial collective modes in 
the BEC-BCS crossover. The data in the upper seven (lower 
eight) rows refer to the sets of measurements taken in the 
shallower (deeper) trap with Uo = 1.8/iK and Ef = 500 nK 
(Uo — 7.3/iK and Ef = 800 nK). The values in parentheses 
indicate la fit uncertainties of individual measurements. Note 
that a systematic scaling uncertainty of ~4% for 1/fcpa results 
from the uncertainty in the atom number N = 2.0(5) x 10 5 . 
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FIG. 3: Normalized compression mode frequency f c versus 
interaction parameter I/Zcfo. The experimental data include 
the small corrections for trap ellipticity and anharmonicity 
and can thus be directly compared to theory in the limit of 
an elongated harmonic trap with cylindrical symmetry. The 
open and filled circles refer to the measurements listed in Ta- 
ble U for lo ± /2tt ps 290 Hz and 590 Hz, respectively. The er- 
ror bars indicate the typical scatter of the data points. The 
filled triangle shows a zero-temperature extrapolation of the 
measurements displayed in Fig. [4] The theory curves refer 
to mean-field BCS theory (lower curve) and QMC calcula- 
tions (upper curve) and correspond to the data presented 
in Ref. [19]]. The horizontal dashed lines indicate the val- 
ues for the BEC limit (f c = 2) and the unitarity limit 
{fc = V 10 / 3 = 1-826). 



BEC regime (1/fcpa > 1) we used the weaker trap with 
loj_/2tt w 290 Hz to minimize unwanted heating by in- 
elastic collisions. Closer to resonance < 1) inelas- 
tic processes are strongly suppressed, but the increasing 
cloud size introduces larger anharmonicity shifts. Here 
we chose the deeper trap with w±/2tt w 590 Hz. On the 
BCS side of the resonance we observed increased dam pin g 
as a precursor of the breakdown of hydrodynamics [9lllll] . 
We thus restricted our measurements to magnetic fields 
below 850 G to ensure low damping rates (j/uj± < 0.01) 
and superfluid hydrodynamics. 

At a given magnetic field, a set of measurements on 
the sloshing and compression modes typically takes a few 
hours. To minimize uncertainties from slow drifts and 
day-to-day variations we always took the sloshing mode 
reference measurement right before or after the compres- 
sion mode data. By repeating measurements under iden- 
tical settings we found a typical remaining fractional un- 
certainty for the normalized compression mode frequen- 
cies of 5 x 10 -3 , which is about 2-3 times larger than the 
fit errors of individual measurements. 

In Fig. [3] we show our final results on the normalized 
compression mode frequency in the BEC-BCS crossover. 
The two theory curves [l9[ correspond to the equation 
of state from mean-field BCS theory (lower curve) and 
the one from quantum Monte-Carlo calculations (upper 
curve). Our data confirm the quantum Monte-Carlo pre- 
dictions and rule out the mean-field BCS theory. In the 



strongly interacting BEC regime (1/fcpa > 1) our data 
are well above the value of 2. This highlights the pres- 
ence of the long-sought beyond-mean-field effects jl7ll in 
collective modes of a strongly interacting gas [3, |26| . 

We finally address the question how non-zero tempera- 
tures influence the compression mode frequency. At uni- 
tarity, a recent experiment [27| has found small frequency 
upshifts with temperature. For a BEC, however, the- 
ory [28j predicts temperature-induced down-shifts, which 
compete with the up-shifts from beyond-mean-field ef- 
fects. We have performed a set of measurements on tem- 
perature shifts in the strongly interacting BEC regime 
(1/fcpa = 0.94). Before exciting the collective oscilla- 
tion, the evaporatively cooled gas was kept in the re- 
compressed trap for a variable hold time of up to 1.5 s. 
During this time residual heating by inelastic processes 
slowly increased the temperature, which we observed as a 
substantial increase of damping with time. The damping 
rate 7 thus serves us as a very sensitive, but uncalibrated 
thermometer @,|27|- Fig.0] where we plot the normalized 
compression mode frequency versus damping rate, clearly 
shows a temperature- induced down-shift. We note that 
previous measurements in the strongly interacting BEC 
regime [tl [ll[ were performed at relatively large damping 
rates in the range between 0.05 and 0.1, where frequency 
down-shifts are significant. 

With our new knowledge on systematic frequency 
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FIG. 4: Normalized compression mode frequency f c versus 
damping rate for 1/kpa = 0.94 ((To = 7.3/TK). The error 
bars represent la fit uncertainties. The dashed lines indicate 
the zero-temperature values predicted by QMC calculations 
(upper line) and mean-field BCS theory (lower line). 



shifts in collective mode measurements, let us comment 
on the previous experiments performed in Innsbruck Q 
and at Duke University [j| [ll[ . We have reanalyzed our 
old data on the radial compression mode and identified a 
previously undetermined cllipticity of e sa 0.2 as the main 
problem in our data interpretation [29| . The fact that we 
had normalized the compression mode frequency to the 
vertical trap frequency (u> c /w y ) led to a substantial down 



shift in the hydrodynamic regime, but not in the colli- 
sionless regime. We furthermore believe that significant 
temperature shifts were present in the previous collective 
mode experiments. In particular for the strongly inter- 
acting BEC regime, temperature shifts in our old data 
on the axial mode [9( and the Duke data on the radial 
mode fTTI j provide a plausible explanation for these mea- 
surements being closer to the predictions of mean-field 
BCS theory than to the more advanced QMC results. 

In conclusion, our work shows that collective modes al- 
low for precision tests of many-body theories in strongly 
interacting quantum gases. In future experiments, the 
observation of collective oscillation modes will serve as a 
powerful tool to investigate strongly interacting superflu- 
ids in a more general context, e.g. in mixtures of fcrmionic 
quantum gases. 
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